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Abstract 

We study asymptotic stability of simple homoclinic cycles in equivariant dynamical 
systems. We introduce a classification of simple homoclinic cycles in M" based on the 
action of the system symmetry group. For systems in M^, we list all classes of simple 
homoclinic cycles. For each class, we derive necessary and sufficient conditions for 
asymptotic stability and fragmentary asymptotic stability in terms of eigenvalues of 
linearisation near the steady state involved in the cycle. 

1 Introduction 

In this paper we continue the investigation ^20j of asymptotic stability of structurally stable 
heteroclinic cycles in a smooth dynamical system 

x = /(x), f:W^W (1) 

equivariant under the action of a non-trivial finite symmetry group F: 

/(7x)=7/(x), forall7eFcO(n). (2) 

Let ^1, . . . £ IR"' be hyperbolic equilibria of ([1]) and Kj : C,j — j- Cj+i) J = 1, ■ ■ ■ ,m, 
(,m+i = ^1, be a set of trajectories from to ^j+i- The union of the equilibria and the 
connecting trajectories is called a heteroclinic cycle. A heteroclinic cycle is structurally 
stable (or robust), if for each j there exists an invariant subspace Pj of ([T]) such that kj 
belongs to this subspace, and in this subspace ^j+i is a sink [21 [131 125]. Cycles where all 
connections are one-dimensional are called simple heteroclinic cycles. 



Heteroclinic cycles often arise in systems related to biology [T71 [TT], fluid dynamics [H 
Uni [21] and game theory [5]. Heteroclinic cycles are often responsible for complex and 
intermittent behaviour [13]. They may have simple geometric structure but complex local 
attraction properties — heteroclinic cycles which are not asymptotically stable can attract 
a positive measure set in its small neighbourhood |3l |5l [TJ [121 [El l22] . In [20] such cycles 
were called fragmentarily asymptotically stable. 

The theory of asymptotic and fragmentary asymptotic stability is yet incomplete even for 
simple heteroclinic cycles. A sufficient condition for asymptotic stability of heteroclinic cycles 
was presented in [H]. Necessary and sufficient conditions for asymptotic stability of simple 
homoclinic and heteroclinic cycles in were proven in [6l[8l[T6|, for fragmentary asymptotic 
stability of simple heteroclinic cycles in in [21] . Asymptotic stability of heteroclinic cycles 
in some particular systems was studied in [S, [5l [6l [71 [HI [121 [151 [221 123] . Necessary and sufficient 
conditions for asymptotic stability of the so-called type A cycles were proven in [H] , and for 
asymptotic and fragmentary asymptotic stability of the so-called type Z cycles in [20] . 

Here, we consider stability of simple homoclinic cycles. We introduce a classification 
of cycles that depends on the action of the system symmetry group. Some our classes 
are comprised of type A homoclinic cycles, some classes of type Z homoclinic cycles, the 
remaining classes of cycles of other types not studied so far. For cycles in we formulate 
necessary and sufficient conditions for stability in terms of eigenvalues of the linearisation 
near the steady state. 

2 Definitions 
2.1 Stability 

Let us recall the definitions of various types of asymptotic stability of an invariant set of 
system ([T|). We denote by $t(x) the trajectory of ([T|) that starts at point x. For a set X 
and a number e > 0, the e-neighbourhood of X is the set 

B,{X) = {x G i?" : rf(x, X) < e}. (3) 

Let X be a compact invariant set of the system ([T]). We denote by Bs{X) its 6 -local basin of 
attraction: 

Bs{X) = {x e : rf(<l>t(x), X) < 6 for any t > and lim rf(<l>t(x), X) = 0}. (4) 

Definition 1 A compact invariant set X is asymptotically stable, if for any 6 > there 
exists e > such that 

B,{X) c BsiX). 

Definition 2 A compact invariant set X is fragmentarily asymptotically stable, if for any 

fxiBsiX)) > 0. 
(Here fi denotes the Lebesgue measure in R".j 
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Definition 3 A compact invariant set X is completely unstable, if there exists S > such 
that fi{BsiX)) = 0. 

We define now various types of stability of a fixed point x of a map (7 : M" — M" (i.e. 
= x). 

Definition 4 A fixed point x G of a map g : — )■ is asymptotically stable, if for 
any 6 > there exists e > 0, such that |y — x| < e implies 

Ig'^y — x| < (5 for all k > and lim g'^y = x. 

Definition 5 ([2UJ, adapted). A fixed point x G of a map g : — )■ is fragmentarily 
asymptotically stable, if for any 6 > the measure of the set 

Vs{x.) =: {y G R" : l^^^ - x| < 5 for all k > and lim g'^y = x} 

is positive. 

Definition 6 A fixed point x G M" of a map g : R" R" is completely unstable, if there 
exists 6 > such that 

Mi^)) = 0. 

2.2 Heteroclinic cycles 

Let ^1, . . . , be hyperbolic equilibria of the system ([T]), ([2]) with stable and unstable mani- 
folds W'^{^j) and W'^{^j), respectively. Assuming ^m+i = (,1, we denote by kj, j = 1, . . . ,m, 
the set of trajectories from to ,^^+1: kj = fl 7^ 0. 

Definition 7 A heteroclinic cycle is a connected component of the group orbit of the equi- 
libria {^1, . . . , C,m} o-nd their connecting orbits . . . , Km} under the action of T . 

The isotropy group of a point x G R" is the subgroup of F satisfying 

Ea; = {7 G r : 7X = x). 

The fixed-point subspace of a subgroup S C F is the subspace 

Fix(S) = {x G R" : ax = X for all a G E}. 

Definition 8 A heteroclinic cycle is structurally stable (or robust), if for any j , 1 < j < m, 
there exist Sj C F and a fixed-point subspace Pj = Fix(Sj) such that 

• ^j+i is a sink in Pj; 

• Kj C Pj. 
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We denote Lj = Pj-i H Pj, by Tj the isotropy subgroup of Lj (evidently, G L^), and 
by the orthogonal complement to Pj in M". 

For a structurally stable heteroclinic cycle, eigenvalues of df{^j) can be divided into four 
sets mUHllIe]: 

• Eigenvalues with the associated eigenvectors in Lj are called radial. 

• Eigenvalues with the associated eigenvectors in Pj^i Lj are called contracting. 

• Eigenvalues with the associated eigenvectors in Pj Lj are called expanding. 

• Eigenvalues that do not belong to any of the three aforementioned groups are called 
transverse. 

Definition 9 (116], adapted). A robust heteroclinic cycle in \ {0} is simple, if for any j 
dim(Pj_i Lj) = 1. 

Definition 10 ^20] A simple robust heteroclinic cycle is of type Z, if for any j 

• dim Pj = dimPj+i; 

• the isotropy subgroup of Pj, T,j, decomposes P^ into one- dimensional isotypic compo- 
nents. 

Note that the first condition in the definition [20] of type Z cycles is redundant, because 
it is implied by the simplicity of the cycles: 

Lemma 1 Let X be a simple robust heteroclinic cycle, i.e. the respective fixed-point sub- 
spaces satisfy dim(Pj_i Lj) = 1 for any j , 1 < j < m. Then dim Pj = dimP^+i for all j , 
1 < j < m — 1. 

Proof: The dimension of the expanding eigenspace for a steady state C,j involved in a 
heteroclinic cycle can not be less than one, and therefore 

dimPj > dimLj + 1 = dimPj_i. 

This inequality applied for each m steady states yields 

dim Pi < dim P2 < • • • < dim P^ < dim Pi 

(recall that the equilibria are cyclically connected, i.e., ^1 = ^m+i)- Here, the leftmost and 
rightmost values coincide, and thus all terms in the inequality are equal. QED 

Since dim(Pj_i Lj) = 1 for a simple heteroclinic cycle, all expanding and contracting 
eigenspaces are one-dimensional. 

Definition 11 A simple robust heteroclinic cycle is of type A', if for any j the isotypic 
decomposition of Pj- under the action of Ej involves only one isotypic component. 

Our A' type cycles are a subset of type A cycles: 
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Definition 12 [T4I IT^ [7^ A simple robust heteroclinic cycle is of type A, if for any j 

• all eigenvectors, associated with X'j, A*, X'j^^ and Xj^i, belong to the same isotypic 
component in the decomposition of Pj- under Sj/ 

• all eigenvectors of df{^j), associated with transverse eigenvalues with positive real parts, 
belong to the same isotypic component in the decomposition of Pj- under Sj. 

Here Xj and A* denote the contracting and transverse eigenvalues ofdf{C,j) with the minimum 
real parts, respectively, and A| the expanding eigenvalue with the maximum real part. 

Note that if a system depends on a parameter and the classification involves conditions 
for eigenvalues, then on variation of the parameter the type of the cycle can change without 
any qualitative change in the overall behaviour of the system^ This is not the case, when 
the classification is based on the action of the symmetry group, as we propose here. 

Definition 13 A heteroclinic cycle is called a homoclinic cycle, if there exists a symmetry 
7 G r such that for any 1 < j < m 

The symmetry 7 is called a twist |2l |25] . 

In this paper we study stability of simple homoclinic cycles, and we use the symbols ^, 
K, P and L without indices. The radial eigenvalues of df{^) are denoted by — r = {—ri}, 
1 < I < Ur, the contracting one by — c, the expanding one by e and the transverse ones by 
t = {ti}, 1 < I < Ut. Here and Ut are the numbers of the radial and transverse eigenvalues, 
respectively. 

3 Classification of simple homoclinic cycles in 

By definition [8], a homoclinic cycle is structurally stable, if there exists a subgroup S C F 
such that the connection from ^ to 7^ belongs to a fixed-point subspace P = Fix(S). We 
introduce a classification of simple homoclinic cycles in M" that is based on the actions of 
S and 7 (note that 7 ^ S) on P^. It is coarser than those proposed in [2^ for homoclinic 
cycles in and employed in [261 E] ioT cycles in and R^. Nevertheless, our classification 
suffices to determine conditions for stability of homoclinic cycles in R^ — for each class, they 
have the form of inequalities for eigenvalues of the linearisation df{C,) (see table 2). 

Let e^, 1 < k < K, he a basis in P^ comprised of eigenvectors of df{^), and let h^, 
1 < k < K, he a basis in P-^ comprised of eigenvectors of dfi^j^). We assume that ei 
is the contracting eigenvector of df{C,), hi is the expanding eigenvector of c?/(7^), and the 
remaining transverse eigenvectors are related: = 'yek, 2 < k < K. Suppose 

P^ = Ui(B...(BUj (5) 

^For instance, cycles of class 3-2-(12)(3) (see subsection 15 . 2 1 and table 2) are not of type A', but for ti < t2 
they are of type A as defined in Jl4^. Such a cycle is asymptotically stable for c < e and ii, ^2 < 0. When ti 
exceeds t2, the cycle just ceases to be of type A, although no bifurcations take place and the stability of the 
cycle does not change. 
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is the isotypic decomposition of P under the action of S. Let the eigenvectors in the basis 
be ordered in such a way that the h^, belonging to one isotypic component, Uj, (whose 
dimension we denote by Ij), have consecutive indices, k = s + 1, . . . , s + Ij; namely, Ui 
is spanned by hi, . . . , h;^, U2 by h/^+i, . . . , h^^+j^, etc. Each isotypic component Uj is also 
spanned by the eigenvectors for some k = iLj^i+i,'iLj^i+2, ■ ■ ■ ,iLj, where we have denoted 
Lj = li + . . . + Ij. We label such a heteroclinic cycle by a sequence of numbers, where 
subsequences associated with individual isotypic components are enclosed in parentheses: 

(zi,i2, • • ■ ,ih){ih+i,ih+2, ■ ■ ■,iL2) ■ ■ ■ {iLj^^+i,iLj^,+2, ■ ■ -^hj)- 

All classes of homoclinic cycles in are listed in table 1. The sequences defined above are 
supplemented by two numbers: the dimension of P"*" and the number of isotypic components 
(e.g. 3-l-(123) labels a cycle with a three-dimensional P"*" comprised of a single isotypic 
component). Note that transverse eigendirections are interchangeable and the resultant 
cycles (e.g., 3-2-(2)(13) and 3-2-(3)(12) ) are identical; thus, only one of them is listed. 



4 Poincare maps for homoclinic cycles in 

Following [T2I [ini [21], in order to examine stability of a homoclinic cycle we consider a 
Poincare map near the cycle. In subsection 12. 21 we have defined radial, contracting, expanding 
and transverse eigenvalues of the hnearisation df{^). Let {u,v,w,z) be coordinates in the 
coordinate system with the origin at ^ and the basis comprised of the associated eigenvectors 
in the following order: radial, contracting, expanding and transverse. 

If 5 is small, in a 5-neighbourhood of ^ the system ([1]) can be approximated by the linear 
systenl§ 

ill = -riui, 1 <l <nr 

^. = (6) 

w = ew 

zi = tizi, 1 <l <nt. 

Here, {u,v,w,z) denote the scaled coordinates {u,v,w,z) = {u,v,w,z)/5. 

^We assume here that all eigenvalues are real. The system under consideration can have a pair of complex 
conjugate radial eigenvalues, if the dimension of the radial eigenspace is larger than one, or it can have a pair 
of complex conjugate transverse eigenvalues, if the homoclinic cycle is of the classes 3-l-(123) or 3-2-(l)(23). 
The radial eigenvalues are not relevant in the study of stability. If transverse eigenvalues are complex, the 
estimates 

fcldzil + \Z2\)\W\-'/' < |zi| < i^idzil + \Z2\)\W\-'/', k2i\z,\ + \z2\)\w\-'/' < \z2\ < K2{\z,\ + \z2\)\w\'*/- , 

where t ~ Re(<i) — Re(i2), can be employed in the proofs of stability and instability, respectively, instead 
of the exact expressions 

zi = zi|w;|-*i/^ Z2=Z2\w\-''/^. 

Since only the values of exponents are important in the proofs, in the conditions for stability (see table 2) ti 
is replaced by Re(ii), and no other modifications are required. 
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Let (uo,fo) be the point in '~f~^P where the trajectory 7~^k intersects with the sphere 
|up + = 1, and q be local coordinates in the hyperplane tangent to the sphere at the 
point (uo,fo). We consider two crossections: 

^(°«*) = {(u,t;,ty,z) : |u|, |z| < l,w = 1} 

and 

^(m) ^ {(ci,w,z) : |q|, \w\, |z| < 1}. 

Near ^, trajectories of system ([T]) can be approximated by a local map (called the first 
return map) : — )■ that associates a point, where a trajectory crosses with 

the point, where the trajectory crossed The global map -ip : — )■ associates 

a point where a trajectory crosses 7if(™) with the point where its previously crossed 
The Poincare map is the superposition g = The w- and z-components of the map g 
are independent of q: this was shown in PH [20] for slightly different systems, but the proof 
can be trivially modified to serve the case considered here. Thus, one can define the map 
g{w, z) that is the restriction of the map g into the {w, z) subspace. The stability properties 
of fixed points of the maps g and g are identical; hence, stability of a cycle is determined by 
stability of the fixed point {w, z) = of the map j?. 

Denote by (j) and ip the restrictions of ip and into P"*". In the leading order, the map </> 

is 

<P{w,{z,}) = ivow'/^,{zs\w\-'^/'}) (7) 
(we use the coordinates {w,z) in and (f , z) in 

Hiout)y ^Q^g ^Y^^^ ^Yie local map IS 
expressed by ([7]) for any homoclinic cycle of whichever class. 

Expressions for global maps are different for different classes of homoclinic cycles. The 
conditions of lemma [2] (see appendix are satisfied for all simple homoclinic cycles in 
except for the 3-2-(l)(23) cycles. We do not consider the 3-2-(l)(23) cycles henceforth in 
this section; the global map for them is derived in appendix |X1 By the lemma, each isotypic 
component of P"*" is an invariant subspace of the map ip. 

Generically, in the leading order the global map is linear in each isotypic component. 
The matrix, C, of the linear map 

(^^«,z^«)=V^(t;«,z«) = c(^^^, (8) 

(here superscripts indicate whether the components are in the basis of eigenvectors of df{C,) 
or of (i/(7^)) is the product C = BA, where A is the matrix of the map ip in the basis of 
eigenvectors of df{^), and B is the matrix of transformation of (f^,z^) into {w'^^,z'^^). In 
the study of stability, we focus on the location of blocks in C = {cij} that vanish because 
the map ip has invariant subspaces. Generically, Cij ^ 0, if and hj belong to the same 
isotypic component in the decomposition Hence, the location of zero entries of C can 
be determined applying the classification presented in section [31 Below we list exhaustively 
the possible forms of matrices C (where non-zero entries are shown by *) for various types 
of homoclinic cycles in M^, and determine the general forms of Poincare maps (see table 1). 
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If the subspace is one-dimensional, then C is an 1 x 1 matrix and the Poincare map 
is just g{w) = Ciiw^/^. 

If the subspace P"*" is two-dimensional, then the classification of simple homoclinic cycles 
in is applicable, since for such cycles is two-dimensional. Alternatively, note that 
such cycles are either of type A' (if the decomposition of P"*" under S has only one isotypic 
component) or Z (if there are two components). In the former case the cycle is classified 
as 2-l-(12), generically none entries of its matrix C vanish, and thus the Poincare map is 
g{w^ z) = {ciiw'^^'^ + C\2z\w\~^l'^ , C2iw'^^^ + C22z\w\~^/^). In the latter case the cycles are of the 
2-2-(l)(2) or 2-2-(2)(l) classes, the matrices of the global map ip are 

(9) 



* \ / * 
* J 1*0 







2-2-(l)(2) 2-2-(2)(l) 

and the Poincare maps are g{w,z) = (cnw^/^, 022-2] w|~*^'^) (for the 2-2-(l)(2) cycle) and 
giw,z) = (ci2z|w|-*/%C2iU7'/') (for the 2-2-(2)(l) cycle). 

If the subspace P^ is three-dimensional, then homoclinic cycles are of types A' (if the 
decomposition (|5]) involves just one isotypic component), Z (if the decomposition involves 
three components), or of other types not studied so far. For a type A' cycle listed in tables 
1 and 2 as 3-l-(123), generically all entries of C are non-zero. For type Z cycles, A is a 
diagonal matrix and P is a permutation matrix (provided vectors e^- and in the bases are 
normalised), and hence C has one of the following forms: 



(10) 

3-3-(l)(2)(3) 3-3-(l)(3)(2) 3-3-(2)(l)(3) 3-3-(2)(3)(l) 

Now suppose decomposition ([5]) involves two isotypic components. Then the cycles are 
either 3-2-(-)(--), or 3-2-(- ■)(■), where dots should be replaced by a permutation of indices 
{1, 2, 3}. The matrices of global maps for such homoclinic cycles have the following forms: 

* 
* 

0**y y*0*y \00*/ \0*0/ \*00/ (11) 

3-2-(l)(23) 3-2-(2)(13) 3-2-(12)(3) 3-2-(13)(2) 3-2-(23)(l) 

Poincare maps for these cycles are listed in table 1. The exponents aj, j =1,2,3, can be 
expressed in terms of eigenvalues of the linearisation df{^) by the relations 

ai = c/e, a2 = —h/e (or —t/e if dimP"*" = 1), = —t2/e. (12) 
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Classification 



Poincare map 



3 
3 

3 
3 



(1)(23) 
;i)(23)m 



2-(2)(13) 
2-(13)(2) 



g{w,zi,Z2 
g{w,u) = 

C22Zl\ 

g{w,Zi,Z2) 
g{w,Zi,Z2) 



w\ 



W'^\ C22Zl\w\''^ + C23Z2\w\'^'-' , C32Zl\w\°-^ + C33Z2 

+ |ti;|«2'=Re((ci2 + ici3)(2i + iz2)^), 



+ C232^l|w|"^ -C23Zi\w\"-'^ + C22Zl\w\"''^ ) 
Ci2Zi\w\°''^ , C2lW°-^ + C23Z2\w\"-^ , C^iW"^ + 0332:2 



W 



Type 



1 - 


1 - 


(1) 




= (CiiW"i) 


A', Z 


2 - 


1 - 


(12) 




z) = (ciiW"! + C122 W °■^ C2lW°-^ + 022^ ^ "^) 


A' 


2 - 


2 - 


(1)(2) 




Z) = {ciiW"'\ C22Z\W\"''^) 


Z 


2 - 


2 - 


(2)(1) 




Z) = {Ci2z\w\"-\ C2iW"i) 


z 


3- 


1 - 


(123) 


C21W 


-21,-22) = (ciiU7"i + Ci2Zl\w\"-^ + Ci3Z2\w\"-^, 

"1 + C22-21 W + C23Z2\w\"-^ , CsiW"-' + C32-2l W + €332:2 W "3) 


A' 


3- 


2 - 


(12)(3) 




Zi,Z2) = (CiiW^i + Ci2-2l W 021^"! + C22-2l W 033^2 ^ ) 





3 


- 2 - 


(23)(1) 


gi 




Zl, 


Z2) 


= (ci2^i w "-^ + 0132:2 


C222:i W "2 + C23-22 W C3ltf"^) 




3 


-3- 


(1)(2)(3) 


g( 




Zl, 


^2) 


= (CiiW"i,C22^l|w|"', 0332:2 


lwl''^) 


z 


3 


-3- 


(1)(3)(2) 


gi 




Zl, 


^2) 


= (CllW;"^C232;2|w^^C322;l 


l^h") 


z 


3 


-3- 


(2)(1)(3) 


gi 




Zl, 


^2) 


= (ci22:i|w|''2,C2iw"S 0332:2 




z 


3 


-3- 


(2)(3)(1) 


gi 


>> 


Zl, 


^2) 


= (cl22:l|^^;|'"^C232;2|^<-'|''^C3l^^7"l) 


z 



Table 1: Poincare maps for different classes of homoclinic cycles in R^. Last column indicates 
cycles of types A' or Z. are the ratios (fT^ of eigenvalues of the linearisation. 

5 Stability 

5.1 Stability of cycles of types A and Z 

In this subsection we review the results on asymptotic stability of cycles of types A [T5] 
and Z [201. 

Theorem 1 (fT^ [73| /. adapted for homoclinic cycles). Let —c, e and tj , 1 < j < J, be the 
contracting, expanding and transverse eigenvalues of df{C,) for the steady states ^ involved in 
a homoclinic cycle of type A. 

(a) If c > e and tj < for all 1 < j < J , then the cycle is asymptotically stable. 

(b) If c < e or tj > for some j then the cycle is completely unstable. 

Stability of the fixed point {w, z) = of the map g associated with a type Z homoclinic 
cycle was studied in [20] by considering the map in the coordinates 

T] = (ln|w|,ln|2:i|,...,ln|2:„J), (13) 

in which the map is linear: 

grj = Mr] + F. (14) 
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Here 





( fli 





. 






(32 


1 


. 


. 


M = B 







1 . 


. 




\ Q-N 





. 


• 1/ 



is the transition matrix of the map i? is a permutation matrix (see section Hj), and the 
entries in the matrix of the local map are 

ai = c/e and a^+i = —tj/e, I < j < J- (16) 

Any permutation is a superposition of cyclic permutations. We assume that vectors ej 
in the basis are ordered in such a way that the first vectors (the subscript s stands for 
significant) are involved in a single cyclic permutation and ei is the contracting eigenvector 
of (i/(0- Matrix B permutes the vectors: ei e2 —i- — ^ ei. Any eigenvalue of 

the upper left Ug x Ug submatrix of M is also an eigenvalue of M, because the upper right 
(J — rig) X Ug submatrix of M vanishes. These eigenvalues are called significant; generically 
they differ from one in absolute value. All other eigenvalues of M are one in absolute value. 
Denote by Amax the largest in absolute value significant eigenvalue of M and by v™'^^ the 
associated eigenvector. Conditions for asymptotic and fragmentary asymptotic stability of 
type Z cycles in terms of Amax and components of v™^^ are stated in theorems [2] and |3l 

Theorem 2 Let M be the transition matrix of a homoclinic cycle of type Z. Suppose all 
transverse eigenvalues of df{^) are negative. 

(a) //|Amax| > Ij then the cycle is asymptotically stable. 

(b) If I Amax I < 1; thcn thc cycle is completely unstable. 

Theorem 3 Let M be the transition matrix of a homoclinic cycle of type Z. The cycle is 
fragmentarily asymptotically stable if and only if the following conditions are satisfied: 

(i) Amax is real; 

Amax > 1; 

(Hi) -ymax^max ^ g ^ I < l,q < N. 

5.2 Stability of homoclinic cycles in 

Conditions for asymptotic stability and fragmentary asymptotic stability for various classes 
of homoclinic cycles are presented in table 2. For type A' cycles the conditions follow from 
theorem [H 
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For type Z cycles the conditions are determined from theorems [2] and [3] by calculating 
eigenvalues and eigenvectors of transition matrices. For the 2-2-(-)(-) cycles the transition 
matrices are 

/ ai \ / (22 1 
y a2 1 J \ ai 

2-2-(l)(2) 2-2-(2)(l) 

The first and second matrices have a one- and two-dimensional significant subspace, re- 
spectively. Calculating the eigenvectors and eigenvalues, we determine the conditions for 
asymptotic stability listed in table 2 (previously found in [161 [201 [21]). When such a cycle is 
not asymptotically stable, it is completely unstable. 
The transition matrices of the 3-3- (•)(•)(•) cycles are 

aiOO\ /aiOO\ /oalOX /aglOX 
02 1 flg 1 ai ag 1 

as 01/ \a2l0/ \a3 1 / \ ai J (18) 

3-3-(l)(2)(3) 3-3-(l)(3)(2) 3-3-(2)(l)(3) 3-3-(2)(3)(l) 

The dimension of the significant subspace of the first, second, third and fourth matrix is 
one, one, two and three, respectively. For the matrices with one- and two-dimensional 
significant subspaces, the conditions of theorems |2] and [3] can be expressed in terms of by 
explicitly calculating the eigenvectors; the respective maps are either asymptotically stable 
or completely unstable. For the fourth matrix, the relations between the entries Oj, that are 
equivalent to the conditions of the theorems, are derived in appendix [Bl Substituting dHj), 
one obtains the conditions listed in table 2. 

Stability of cycles which are not of types A' or Z was not yet studied. In appendix O 
we derive conditions for asymptotic and fragmentary asymptotic stability of the 3-2-(2)(13) 
cycle. Its Poincare map is 

^(^,^1,2:2) = (cl22;l|w|"^C2lw"l + C23Z2\w\''' , c^iw"^ + €332:2 |w; 1"^) . 

If w, Zi and Z2 are small, then in the domains w°'^~"'''^ ^ Z2 and w°'^~"'^ ^ Z2 (which we 
denote by flj and fljj, respectively), the map g is approximated by the maps 

gi{w,zi,Z2) = {ci2Zi\w\"-\ C2iw''\ csiw"'^) and gn{w,zi,Z2) = (ci22;i|w|''^ C23^2k^^ 033^21^!"^), 

respectively. Let M/ and Mu denote the transition matrices of the two maps, 

/ as 1 \ / a2 1 \ 

= ai , M,, = a3 1 , 
\aiOO/ \a3 01/ 

and A"^"^ their largest in absolute value eigenvalues, and v^^ the respective associated 
eigenvectors, 

= (1, A^ - as, A^ - as) and v^^ = (1, A^^^ - a2, A^-^ - a2). 



n 7\ 
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and A2 3 and \2^.^ the remaining eigenvalues. 

Theorems S] and O (see appendix E]) can be summarised as follows: 

(i) the cycle is asymptotically stable if and only if all entries of the transition matrices M/ 
and Mil are non-negative and min(|A'^|, |A^^|) > 1; 

(ii) the cycle is fragmentarily asymptotically stable if and only if at least one of the following 
conditions is satisfied: 

- the matrix M/ satisfies the conditions of theorem [3] (namely, is real, A''^ > 1, 
A^ > I min(A2, Ag)] and vlvl > for all 1 < fc, / < 3) and A^ - 02 > (ai - a^)'^; 

- the matrix Mjj satisfies the conditions of theorem [3] and A^^ — 02 < (ai — 03)^^; 

(iii) if the cycle is not fragmentarily asymptotically stable, then it is completely unstable. 
We use the T;-coordinates (fT3l) to define the set 

= {('7i,'72,?73) : 773/771 > ai - as + e}. 

For any e > the set flj^^ is a subset of (upon the change of coordinates ( fT3l) ). Note 
that the condition \^ — a2 > (cti — 03)"^ for fragmentary asymptotic stability (see (ii)) is 
equivalent to the condition that the eigenvector belongs to Qj^^ for some e > 0. An 
analogous equivalence holds for the matrix Mu and the set 

^ii,e = {iVi,V2,V3) ■ rj^/rji < ai- a^- e}. 

Conditions (i) and (ii) expressed in terms of eigenvalues of df{E,) are listed in table 2. For 
other cycles from the 3-2-... classes, the results can be obtained similarly — by introducing 
analogous matrices M/ and M// and applying to them the results of theorems [2] and [3] 
together with the conditions that C and v'^'^ C for some e > 0. 

6 Conclusion 

We have introduced a classification of simple homoclinic cycles in and identified all classes 
of homoclinic cycles in M^. For each class in we have derived necessary and sufficient 
conditions for asymptotic and fragmentary asymptotic stability in terms of eigenvalues of 
linearisation near the steady state involved in the cycle. Transition matrices were used in 
the derivation of conditions for stability for cycles from the new classes, as in [20] for type Z 
cycles, although now the conditions for stability and the derivations are more involved. The 
approach based on the use of transition matrices is likely to be fruitful for investigation of 
stability of other simple cycles, i.e., for homoclinic cycles in for n > 5 and for heteroclinic 
cycles in M" for n > 4. 

In M^, our classification is sufficient to determine the stability of a homoclinic cycle: the 
stability is controlled, in a class-dependent way, by eigenvalues of the linearisation. It is 
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Classification Conditions for stability 



1 


- 1 - 


(1) 


A. 


s.: 


c > e 




2 


- 1 - 


(12) 


A. 


s.: 


c> e, t < 




2 


- 2 - 


(1)(2) 


A. 


s.: 


c> e, t < 




2 


- 2 - 


(2)(1) 


A. 


s.: 


c-t>e, t <0 




3 


- 1 - 


(123) 


A. 


s.: 


c> e, ti < 0, t2 < 




3 


- 2 - 


(12)(3) 


A. 


s.: 


c> e, ti < 0, ^2 < 




3 


- 2 - 


(1)(23), 


A. 


s.: 


c> e, ti < 0, ^2 < 




3 


- 2 - 


(l)(23)m 


A. 


s.: 


c> e, ti < 0, t2 < 




3 


- 2 - 


(2)(13) 


A. 


s.: 


c - ti > e, ti < 0, t2 < 










F. 


a. 


s.: ti(c - ti) + (c - t2)^ - ce > 0, ti 


< e, t2 < 


3 


- 2 - 


(13)(2) 


A. 


s.: 


c> e, ti < 0, t2 < 










F 


. a. 


s.: c(c — ^2) — e^ — tie < 0, ti < e 
or c(c — t2) — e^ — tie > 0, ti 


, t2 < 0, ti + t2 < 

< e, c > e 


3 


- 2 - 


(23)(1) 


A. 


s.: 


c — max(ti, t2) > e, ti < 0, t2 < 




3 


-3- 


(1)(2)(3) 


A. 


s.: 


c> e, ti < 0, t2 < 




3 


-3- 


(1)(3)(2) 


A. 


s.: 


c> e, ti < 0, t2 < 




3 


-3- 


(2)(1)(3) 


A. 


s.: 


c - ti > e, ti < 0, t2 < 




3 


-3- 


(2)(3)(1) 


A. 


s.: 


c - ti - t2 > e, ti < 0, t2 < 










F. 


a. 


s.: c — ti — t2 > e, tit2 + ce > 0, ctf 


+ et| < 



Table 2: Conditions for asymptotic stability and fragmentary asymptotic stability of different 
classes of homoclinic cycles in in terms of eigenvalues of the linearisation df{C,). 

unclear whether the classification suffices to determine conditions for stability in higher- 
dimensional systems. Note that a finer classification involving structure angles was proposed 
for homoclinic cycles in and [25l [26] . 

A study of stability is often followed by a study of bifurcations. A natural continuation 
of the present paper is an investigation of resonance bifurcations of simple homoclinic cycles 
in M^, similar to the investigation [7] for homoclinic cycles in R^. 
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A Global map for the 3-2-(l)(23) homoclinic cycles 

Consider a map : — )■ M'^ equivariant under a symmetry group S C 0(3). Suppose 

(a) S decomposes M.^ into two isotypic components 

where the dimension of f/i is one and of f/2 is two; 

(b) for any x G there exists cr G S such that crx 7^ x. 
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[16 

[17: 

[18; 
[19; 

[20; 

[21 
[22; 
[23; 
[24 

[25; 
[26" 



Consider the {x,z) coordinates in M^, where x is the coordinate in Ui and z in 1/2- In this 
appendix we determine the leading terms of the expansion of ip in small x and z. We use 
the following lemma |2U] : 

Lemma 2 Let a group S act on a linear space V . Consider the isotypic decomposition of 
V under the action ofL: 

V = Uo®Ui® ...®Uk. 

Suppose 

• the action of S on Uq is trivial; 

• any cr G S acts on a Uk, ^ ^ k < K , either as I or as —I. 

Then for any collection of subscripts 1 < ii, ■ ■ ■ ,ii < K there exists a subgroup Gi^^,,,/i_^ C S 
such that the subspace 

is a fixed point subspace of the group Gij^^__,^i^. 

By (a) and (b), there exists a symmetry ai G S such that cri(a;,0) = (— x, 0). S can act 
on U2 in three ways: 

(i) there exists (J2 G S such that (T2 7^ erf for any k and o"2(0, z) = (0, —z); 

(ii) no (72 satisfying (i) exists; there exists as G S such that 0-3(0, z) = (0, e^'^'/'^z), where 
k > 1 is odd, and ai{x, z) = (— x, z); 

(iii) no (72 satisfying (i) exists; there exist (T3 G E such that (73(0,2;) = (0, e^'^'/'^^;), where 
A; > 1 is odd, and <Ji{x, z) = (— x, —z). 

Lemma 1 implies that in case (i) the x— and 2;-components of are 

ip'^ = xF{x, z, z), = zG{x, z, z) + z^H{x, z, z), 

where F is real, generically F(0, 0, 0) 7^ and G{0, 0, 0) 7^ 0, s > is odd (it is determined 
by S). In cases (ii) and (iii) the components of ip can be calculated by simple algebra: 

(ii) ijj^ = xF{x, z, z), if)^ = zG{x,z,z) + z^^^H{x,z,z) 

(iii) = xF{x, z,z) + z^ J{x, z,z) + z^ J{x, z,z), ip^ = zG{x, z, z) 

where F is real, generically F(0, 0, 0) ^ 0, G(0, 0, 0) ^ and J(0, 0, 0) ^ 0. Thus in cases 
(i) and (ii), for small x and z the asymptotically largest terms of ip sxe 

ip'' = ax, %p' = hz + cz% (19) 

where a is real, h and c are complex, and generically a 7^ and 6 7^ 0. (Here h and c 
are further restricted by the action of other symmetries from S, but these restrictions are 
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insignificant for tlie study of stability of the Poincare map.) In case (iii) the asymptotically 
largest terms of the map ip are 

= ax + bz'' + bz'', iP' = cz, (20) 

where a is real, b and c are complex, and generically a ^ 0, b ^ and c ^ 0. We refer to a 
cycle with the global map (|T9ll in cases (i) or (ii) as a 3-2-(2)(13) cycle, and to a cycle with 
the global map ( 120|) in case (iii) as a 3-2-(2)(13)m cycle. 

B Stability of the 3-3-(2)(l)(3) homoclinic cycle 

In this appendix we derive necessary and sufficient conditions for asymptotic stability and 
fragmentary asymptotic stability of the 3-3-(2)(l)(3) homoclinic cycle in terms of eigenvalues 
of the linearisation near the equilibrium. As noted in section |U such a cycle is of type Z. 
Conditions for stability of type Z cycles in terms of eigenvalues and eigenvectors of their 
transition matrices are given by theorems [2] and [3l The transition matrix of the cycle is 

/ a2 1 \ 

M= as 1 (21) 
\ai / 

(see f|T8l) ). where Oj are related to eigenvalues of the linearisation by f|T6|) . ai > (a2 and as 
can have arbitrary signs). 

Let Ai, A2 and As be the eigenvalues of ( 12T1) : Ai denotes the largest eigenvalue if all 
eigenvalues are real, or the real eigenvalue if the matrix has complex ones. Let w denote 
the eigenvector associated with Ai. By theorem |2l the necessary and sufficient conditions for 
asymptotic stability are 

a2 > 0, as > 0, max \Xj\ > 1. (22) 
j 

By theorem [31 the cycle is fragmentarily asymptotically stable if and only if 

Ai > 1; (23) 

Ai > max(|A2|,|As|); (24) 

WiWj > for any 1 < i, j < 3. (25) 

By applying the following lemmas, one can avoid calculating the eigenvalues Aj by Car- 
dano's formulae for the roots of a cubic polynomial. 

Lemma 3 Let all > in matrix l[21\) . Then maxj |Aj| > 1 if and only if 

ai + a2 + as > 1. (26) 
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Proof: Eigenvalues of matrix f l2T|) are roots of its characteristic polynomial 

Pm(A) = -A^ + aaA^ + asA + ai- (27) 

Suppose the inequality fl26|) is satisfied. This implies Pm(1) > 0. Since pm{oo) < 0, the 
polynomial pm has a root larger than one, and thus maxj |Aj| > 1. 

We prove now the converse. Denote by Amax the maximal in absolute value root of Pa/(A). 
Since all a,- > and lAmaxI > 1, 



as ai 
ai + a2 + da > 02 + 7-; ; + 7-; ^ > 



I Amax I I Amax 



02 + -; h 



/^max ^max 



Amax I ^ 1 • 



QED 

Components of the eigenvector w associated with the eigenvalue Ai satisfy the equations 

02^1 + ^2 = AiWi, (28) 
a^wi + W3 = X1W2, (29) 
aiWi = X1W3. (30) 

By the Viete formulae for the roots of the characteristic polynomial pm, 

A1 + A2 + A3 = 02, (31) 
— A1A2 — A2A3 — A1A3 = 03, (32) 
A1A2A3 = ai. (33) 

Lemma 4 Eigenvalues and eigenvectors of matrix f21\} satisfy conditions ( fl^) -( fl3]) if and 
only if the following four inequalities hold true: 



ai > 0, 

ai + a2 + 03 > 1, 
0203 + ai > 0, 
aial + a| > 0. 



(34) 



Proof: The cubic polynomial pm has either three real roots or one real root and two complex 
conjugate ones. We consider the two cases separately. 

Suppose all eigenvalues Aj are real. 

Assume that (I23l)-(l25l) hold true. By virtue of ( l30l) and ( l23i) . and since wi and W3 have 
same signs ( 125|) . we have ai > 0; hence by ( l33l) A2 and A3 have same signs. Equations (l3Ti) 
and ([28]) yield 

- Wi{\2 + A3) = W2, (35) 

whereby A2 + A3 < (see fl25]) ). Therefore, A2 < and A3 < 0. Thus, and the identity 

(A? - A2A3)(A2 - A3Ai)(A^ - A1A2) = a^al + 4 (36) 
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(which follows from the Viete formulae) yield 

aial + 03 > 0. 

The characteristic polynomial pm (1271) has only one positive root that is larger than one. 
Since pm{oo) < 0, this implies Pm(1) = — 1 + cti + 02 + as > 0. The identity 

(Al + A2)(A2 + A3)(A3 + Al) = -0203 - ai, (37) 

and yield 0203 + Oi > 0. Thus, all inequalities in ( IMl) are proven. 

We prove now the converse assuming that (15^ is satisfied. Since ai > 0, by virtue of (15^ 
Al > for our ordering of Aj, and A2A3 > 0. The inequality 0203 + ai > and fl37|) imply 
that A2 < and A3 < 0. From fl30l) we deduce that wi and W3 have same signs; by virtue of 
fl35|) Wi and W2 also have same signs, which proves fl25|) . Since PAf(l) > and pa/(oo) < 0, 
fl25]) holds true. Due to f l5B]l and f p7|) . the inequalities 0203 + ai > and aia^ + al > imply 
that the condition fl2^ holds true. 

Suppose now the polynomial pm has one real root Ai and two complex conjugate roots 
A2,3 = a ±if3. Identities (l37|l and (136|) yield, respectively, 

0203 + ai = -2a((Ai + a)^ + /3^) (38) 

and 

aia^ + a^ = e(A?-A2A3)(a2 + /32)-^ (39) 

where 

e = (Ai(a^ + (3^) -a^ + 3al3^y + {3a^(3 - {3^ > 0, (40) 

unless Oi + 02 + 03 = and Sa"^ = 

By fl39l) and fHOl) . condition flM|) is equivalent to the inequality aia| + ag > 0. Since 
A2A3 = IA2P, fl30|) and fl33l) imply that wi and W3 have same signs. By fl38l) the inequalities 
0203 + ai > and a < are equivalent. By virtue of fl3Tl) . fl28|) reduces to 

— 2a;u;i = 1^2, 

and hence 0203 + ai > is equivalent to (125|) . 

For the characteristic polynomial ( 1271) with only one real eigenvalue the conditions 
( l23l) and 01 + 02 + 03 > 1 are equivalent. Finally, ( l33l) imply that Oi and Ai have same signs. 
QED 

Substituting expressions (lT6l) into the inequalities in the statements of lemmas 13] and HI 
we establish conditions for stability of the cycle presented in table 2. 
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C stability of the 3-2- (2) (13) homoclinic cycle 

In this appendix we derive necessary and sufficient conditions for asymptotic and fragmentary 
asymptotic stability of the 3-2-(2)(13) homochnic cycle. The Poincare map near the cycle is 

g{w,zi,Z2) = iAzi\w\''\ Bw"' + C Z2\w\''\ Dw"' + Ez2\w\''^) (41) 

(see table 1). Recall that ai = — c/e, and therefore by the definition of expanding and 
contracting eigenvalues (see subsection 12.21) oi is always positive. 

If w, z\ and Z2 are small and Z2 ^ then the map ( HTl) can by approximated by 

gi{w,z^,Z2) = {z^\wY\w''\w'''). (42) 

If w, Z\ and Z2 are small and Z2 ^ w°'^~°''^ (which is possible only for ax > a^), the map can 
be approximated by 

gn{w, zu Z2) = {zx\w\''\ Z2\w\''\ ^2^1"'). (43) 
Consider the maps gi and gn restricted to In the new coordinates 

T] = (In w, In , In 2:2) , (44) 

the maps gi and gn are linear. The matrices of the two maps, 



M/ = ai (45) 



and 





M,, = aa 1 , (46) 



have eigenvalues 



and 



A{,2 = ^(«2 ± {4 + 4ai)i/2) and = 0, (47) 



1 



Mi = + «2 ± ((1 + a2f + 4(03 - a2)Y'^) and A^^ = 0, (48) 
respectively. Since ai > 0, Af 2 fHTl) are real. Eigenvalues Af^2 real if and only if 

(1-02)^ + 403 > 0. 

We assume now that (1 — 02)^+403 > and denote the largest eigenvalues of the matrices 
Ml and Mu by A''^ and A^^. The associated eigenvectors are 

= (l,A^-02,A^-02) and v^^ = (l,A^^-02,A^^- 02), (49) 



20 



respectively. Simple algebra reveals that the inequality 



is equivalent to 

02(01 — 03) + (oi — 03)^ — Oi < 0. (50) 

It is straightforward to establish that ( l50i) is not satisfied, if (1 — 02)^ + 4a3 < 0. 

We present here proofs of theorems and lemmas only when fl50l) is satisfied. (When 
the converse holds true, the proofs remain similar.) For convenience, without any loss of 
generality we also assume 

A, C, E > 0, B,D<0 and - C/B > -E/D. (51) 

There exists p, < p < 1, such that for any and p" from the interval [— eigenvalues 
of the matrices 

a2 l+p' 0\ / a2 l+p" 

Mi^p, = \ ai and Mn^p. = 03 1 | , (52) 
ai / ' \ 03 1 

are real, and the largest eigenvalues (which we denote by X^, and X^i, respectively) satisfy 
the inequahties 

Aj, < A^^, (1 - A0(1 - A^,) > and (1 - X^^){1 - X^i) > 0. (53) 
Let s > satisfy 

s<miJ ^'^^~P\ Y (54) 

V fli - 03 ' {ai - %)(! -p)J 

For ai > 03, we partition into five regions (if ai < 03, such partitioning is not needed, 
see the proof of theorem H]): 

ill = {{w,z,,Z2):\w\^'-'''>-^\z2\{l + \z2n} 

Jd 

flu = {(u;,zi,z2) : < -^1^21(1- k2n} 

^III = {{w,ZuZ2):-^Z2{l-\z2n<W^^~^'<-^Z2{l + \z2n} 
fllV = {{w,Z,,Z2):-^Z2{l-\z2n<W^^-'''<-^Z2{l+\z2n} 

D B 

VLv = {{w,ZuZ2) : - — \Z2\{1-\Z2\') KW^'-"' < - — Z2{l+\Z2\'), {w , Zi, Z2) ^ ^jjj U ^jy} . 

(55) 

Define the sets 

u = {{w,Zi,Z2) : |z2r+^ < \zi\ < \Z2\'-''}, (56) 
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and 

= w n (^j, 1<J<V. (57) 
(jJi and (jJii are partitioned further: ui = ujij U Uiji and uju = ^ii,i U Uu^i, where 

uj J = n g~^0Jj and uj jj = uj\ uj j] 

(58) 

uJiiji = uJii n fi-^^w// and uiij = un \ uuji. 

Denote 

Q(/i, /2, /i) = {{w, zi, Z2) : fi{z2) <w < 72(^2), 2:1 = h{w, Z2)}, (59) 

^) = ^2) : w = f{z2), zi = h{z2)}, (60) 

by 5e the e-neighbourhood of the point {w, zi, Z2) = and 

fij(e) = n 5„ coj{e) =cojn B,, Qif,, /a, h; e) = g(/i, /a, /i) n B,, q{f, h; e) = q{f, h) n 5,. 
Furthermore, 

/i(x) ~ /2(a^) denotes that /i(x) — /2(a;) = o(/i(x)) when x — > 

and 

/i(x) ~ f2{,x) denotes that /i(x) ~ Ff2{x) for a constant F 7^ 0. 

For ai > max(0,a3), the following properties of the map g are direct consequences of 
6ID: 

(a) g is one-to-one, except for the hyperplane w = 0; 

(b) if q{f,h) C Uj and / ~ Z2, then gq{f,h) = q{f,h) where / ~ Z2, /i ~ ^2 and 
(l — p + a2a)ai^a^^ < a < (1 + p + a2a)a]^^a''^; 

(c) if q{f,h) C OJjj and / ~ Z2, then gq{f,h) = q{f,h) where / ~ z^, h Z2 and 
{1 — p + a2(y){l + a^a)^^ < a < {1 + p + a2«)(l + 03«)~^ 

(d) if q{f, h) C uin and / - (^5-1^2)^/^"'""'^ ~ ^2 (note that a > s + (ai - 03)-^), then 
5'?(/, h) = q{f, h), where / ~ z^, /i ~ z^, /3 = (ai - 1 + a(ai - a3))ar^ 

and ((1 - p){ai - 03) 02)0^"^ < a < ((1 +p)(ai - 03) + 02)0^"^; 
g'^lif, h) = q{f, h) where / ^ z!^,hr^ Z2 and 

ai — 1 -|- a(ai — 03) 02 Oi — 1 + a(ai — 03) 02 

H < a < — — : + 



ai((l +p)(ai - 03) + 02) ai ai((l - p)(ai - 03) + 02) cti 
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(e) if q{f, h) C w/y^and / - (ED-^^aJ^/^"'""'^ ~ (note that a > s + (ai - ag)"^), then 
5'Q'(/) ^) = Q'(/) h) where / ~ h Z2, P = ai{ai — 1 + a{ai — 03))"^ and 

{{l-p){ai-a3)+a2){ai-l+a{ai-a3)y^ < a < {{l+p){ai-a3)+a2){ai-l+a{ai-as)y^; 

g'^q{f, h) = q{f, h), where / ~ ; ^ ~ -^2 and 

1 flo 1 Oo 

— <a< —— ^ + — ; 



ai((l +p)(ai - as) + 02) ai ai((l - p)(ai - 03) + 02) ai 

(f) if q{f, h) C cjy, then gq{f, h) = q{f, h) where / ~ h Z2, 1— p<P<l+p and 



-1. 



((1 - p){ai - as) + a2)a^ ^ < a < ((1 + p){ai - as) + 02)0^ 

(j) if QUij2.h) C ui, fj ~ 4 and /i - /2 ~ F^zl then gQ{h,f2,h) = Q{fi,f2,h), 
where /i — /2 ~ ci2F fizl and 7 = 7a]'"'^a~^; 

(i) if Q{fi,j2, h) C w//, /i ~ z^ and /i - /2 ^ i^/i^L then c/Q(/i, /2, /i) = Q{fi, h), 
where /i — /2 ~ CL2FfiZ2 and 7 = 7(01 + aso)"^] 



By definitions of the sets fij, for s satisfying ([5i 

51(^7 U fi// U fiy) Cw. (61) 



If a2 > and fl50|) holds true, then the hsted above properties of the map g imply that, for 
a sufficiently small e, the following inclusions take place: 



g{fliii{e)) n g{Qiv) C w/, g^Qivie)) fl 5'(^^/y) C w/, 
g{ujij{e)) C w/,/, fi'(wy(e)) C W/,/, ^(^(^/^(e)) C W/,/; 
if 02 < and (l50l) holds true, then, for a sufficiently small e, 

g{njjj{e)) n fi'(^^m) C w//, g{nivie)) n 5'(^^m) C W//; 
g{uJii,ii{e)) C w//,//, 5((a;y(e)) C W//,//, g^{uiii{e)) C c^//,//. 



(62) 



(63) 



For a set X C M'^ we denote by ^zi{X) the measure (in M^) of the projection of the set 
X on the plane Zi = {{w, Zi, Z2) being cartesian coordinates in R'^). It is easy to show that 
for fl ~ Z2, /i — /2 ~ -^"-22! 7 > « and a small e we have 

r (7 + l)-iFe^+i ,a > 1, 

/^..(Q(/i,/2,/^;e))^ ^ ^ (64) 

[ (7 + l)-iFe(^+i)/" ,a < 1. 
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Lemma 5 Let Qq, Qi, . . . , Qk be the following sets: 

Qk = Qifki, fk2, h), fki ~ fk2 ^ z^, 0<k<K, foi - /o2 ~ Bz^, 

Qk C Qo, fki - fk2 ~ Cfc^^^ l<k<K, 
where P > a, 'yk> f3, Ck = 7fc//3 and (3 > /3m- 
For any d > there exists 

e = e{B, Ci,..., Ck, /3m, ci, . . . , Ck, d) > 

such that 

K 

J2 i^zr (Qk nB,)< dfi,, (Qo n 5,) . 

k=l 

Proof: We present here a proof for a > 1 (the proof for a < 1 is similar). By ([& 

/i,, (Qo n B,) ^ B{^k + 1) B{ck(3^ + 1) 

Any e, < e < 1, such that 

. /rf g(c,/3m+i) y^^^-^^^'^^^ 

satisfies all conditions of the lemma. QED 
Lemma 6 Let the map r : M — )■ M, 

, , ax + b 

= 

cx + d 

be defined for all x G [xi,X3], where X3 > Xi. Denote the superposition of k (identical) maps 
r by 

r^{x) = r o . ^ . o r (x) . 

k times 

(a) If r{xi) > Xi, r{x3) < X3 and r{x2) < X2 for some X2, Xi < X2 < X3, then 

X — r{x) > min(x2 — r{x2),X3 — r(x3)) for any x G [x2, X3]. 
Hence, there exists K > such that 

r'^{x) < X2 for any x G [0:2, X3] and for some < k = k{x) < K. 

(b) If r{xi) < Xi, r{xz) > X3 and r{x2) > X2 for some X2, Xi < X2 < Xs, then 

r{x) — X > min(r(x2) — X2, t(x3) — X3) for any x G [x2, x^]. 
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Hence, there exists K > such that 

r^{x) > X3 for any x E [x2, X3] and for some < k = k(x) < K. 

(c) If r{xi) < Xi and r(x3) > X3, then there exists X2, Xi < X2 < x^, such that r{x2) = X2 
and 

for any x G [xi, X2[ there exists k = k{x) > such that r'^{x) < xi, 
for any x e]x2,X3] there exists k = k{x) > such that r^{x) > X3. 

(d) If r[xi) > Xi and r^x^) < X3, then 

r{x) G [xi, Xs] for any x G [xi, X3]; 
there exists X2, Xi < X2 < X3, such that 

r(x2) = X2 and lim r^[x) = X2 for any x G [xi,X3]. 

fc— >oo 

Proof: We present here a proof for ) (the proofs for other cases are similar). The 

function 

, , , , . ah — adc^^ 

n{x) = X — r[xj = X — 

c cx + d 

has either no local extrema, or two local extrema Xextr satisfying 

(cXextr + d)'^ + cb — ad = 0. 

In the former case, h is a growing function on each interval of its domain, and the statement 
is trivial. In the latter case, the local extrema are separated by the point x = —dc~^, 
where r(x) is undefined. Hence, any interval, where h{x) is defined, contains at most one 
local extremum. If a local extremum belongs to [x2,X3], the extremum is a local maximum 
because of the condition h{x2) > h{xi). Therefore, h{x) takes the minimum over [x2,X3] at 
an endpoint of the interval. 
Set 

C = — min(r (X2) — X2, r (X3) — X3) > and K = [(X3 — X2) /C] + 1. 
(Here [q] denotes the largest integer not exceeding a real number q.) For any x G [x2,X3] 

X2 — r^(x) = (x2 — x) + (x — r(x)) + . . . + {r^^^{x) — r^{x)) > (x2 — x) + KC = X3 — x > 0. 

QED 

In the following two lemmas we reformulate (omitting proofs) the results of pO] on 
asymptotic stability of the fixed point x = of a map g : M.-^ ^ M.'^ , 

g{x) = (|xi|''"|x2ri^ ■ ■ ■ |xj|"", . . . , \xi\'''^X2\'''^ ■ ■ ■ IxjI"'') (65) 
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in terms of eigenvalues and eigenvectors of its transition matrix 

(an ... aij \ 
. (66) 
aji ... ajj J 

To prove stability, the coordinates 

(?7i, ...,r]j) = (In . . . ,ln \xj\) 

were introduced in [20]. In the new coordinates the map g becomes linear, ^(x) = Mr]. The 
initial condition r] was decomposed in the basis of eigenvectors of M, denoted by Cj-: 

J 

ri = Y.hCr (67) 
i=i 

Let Amax and C,^^^ denote the largest in absolute value eigenvalue of M and the associ- 
ated eigenvector, respectively. The proof exploited the fact that for A; — )■ oo the iterates 
M^r) become aligned with ^max^^axC™'^'') where 6max is the factor in front of C,^^^ in the 
decomposition (l67|l . 



Lemma 7 Let all eigenvalues of the transition matrix ( loq) of the map g ( Ip3]) differ from 
one in absolute value and all entries he non-negative. 

(a) If |Amax| > 1; then the fixed point x = of is asymptotically stable, i.e., for any 
(5 > there exists e > such that any x, |x| < e, satisfies 

15^'' (x) I < 5 for all k>0 and lim g''(x.) = 0. 

(b) If I Amax I < 1; then the fixed point x = ofgis completely unstable, i.e., there exists 
6 > such that for any x 7^ there exists k > satisfying 

\g'{^)\>6. 

Lemma 8 Let the absolute values of all eigenvalues of the transition matrix of the map g 
/ fg3]) differ from one. 



(a) If 

Amax IS real, A^ax > 1 and C^'^^^CT^ > for alll<l,s< J, (68) 

then the fixed point x = 0/ the map g is fragmentarily asymptotically stable, i.e., for 
any 6 > the measure of the set of points x such that 

|5^^'(x)| < S for allk>0 and lim ^^^(x) = (69) 

is positive. For a sufficiently large C > 0, the condition is satisfied by all points 
in a small neighbourhood of —CC,^^^ . 
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(b) If at least one of the conditions in (E^i ^■^ iT'Ot satisfied, then the fixed point x = 0/ 
the map g is completely unstable, i.e., there exist 6 > and, for almost all x (for all x 
except for a set of a zero measure), k = k{x) > such that 

\g\^)\>5. 

Lemma 9 Suppose the inequalities 

> 1, A^^ > 1, ai > 0, as > 02 > 

hold true for matrices Mj ( f73] j and Mjj ( f^^ . (The first inequality implies that 01 + 02 > 1.) 
Then the largest eigenvalue A^'^^ of the product MjMjj is real and 

X'''' > 1. 

Proof: The characteristic equation for the matrix MiMu is 

A'^ — A^(a2 + 03 + ai) — Xai^a^ — 02) = 0. 

One root is zero and the remaining two satisfy 

A1A2 = 0-1(0.2 ~ and Ai + A2 = O2 + 03 + Oi. 

By conditions of the lemma, the product is negative, and the sum is positive and larger than 
one. Therefore, Ai and A2 are real, A2 < and Ai > max(l, IA2I). QED 

We say that x < y, where x, y G M'^, if Xj < yj for any 1 < j < J. 

Definition 14 A map g : — > M"' is bounded from above (from below) in V G M"^ by a 
map : R-^ — > R-^, ifdi'^) > fi'(x) (^(x) < g{'x.), respectively) for any x G V. 

Lemma 10 Suppose all entries of the transition matrix M of the map g ( fg3]) are non- 
negative. If \g\ is bounded from above (below) by g in V, then \g^\ is bounded from above 
(below) by g^ in W G V , where 

W = {xeV : /(x) G V and |^*^(x)| E V for all < k < K}. 

Proof: We present here a proof for \g\ bounded from above by g (the proof for \g\ bounded 
from below by g is similar). Note that if |x| < |y|, then 

|^(x)| <5^(x) <5^(y). 

We prove that 

\g''i^)\<n^) (70) 
by induction in K. Suppose fITOl) holds true foi I = K — 1. Then 

\gg^-\^)\<g\g^-\^)\<n^). 

QED 
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Theorem 4 Consider the map 

g{w,Zi,Z2) = (Mkr',5u7"i +C^2k|"^L'w"^ + Ez2k|°^), (71) 
where ai > and the exponents ai, 02 and 03 satisfy the inequality (5^) . 

(a) If 

02 > 0, as > and ai + 02 > 1, (72) 
then the fixed point {w, zi, Z2) = of the map g is asymptotically stable. 

(b) If 

02 < or 03 < or ai + 02 < 1, 
then the fixed point {w,zi,Z2) = of the map g is completely unstable. 

Proof: Define the maps 



gd,i{w, Zi, Z2) = {\zi\\w\"-'^'^''-,\w\"''^'^'^,\w\"-'^'^^) and 



(73) 



whose transition matrices are, respectively, 

a2 + d 1 \ / a2 + d 1 

Mdj = I ai + and M^ji = + d 1 | . (74) 
ai + d / \ as + d 1 

If Oi < 03, then for smaU w, Zi and Z2 the map ( 17Ti) can be approximated by gi (i.e., for 
a sufficiently small e and arbitrarily small \d\, the map g is bounded in from above by 
gdj for (i < and bounded from below by gdj for c? > 0). The theorem follows from lemma 
[TOl In what follows, we present the proof for the case ai > 03. 

(a) Recall that (l50l) and (!72l) imply that A77 is real and Xu > A/ > 1. Also recall that a 
fixed point (w, zi, ^2) = of a map g is asymptotically stable, if for any 5 > there exists 
e > such that 

\g^{w, zi, Z2)\ < 5 for all A; > and lim g''{w, z\, Z2) = (75) 

for any {w, Zi, Z2) G 5^. 

We show that the iterates g^{w^ Z\, Z2) are bounded from above by constructing different 
bounds for the cases 02 < 03 and 02 > 03. If 02 < 03, we show that 

\g\w, zi, Z2)\ < ga^'iiigdjgdJif'gd'iiw, zi, Z2) where fci + 2k2 + h = k (76) 
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for any {w,Zi,Z2) G -Be(d), any d < and a sufficiently small e{d) > 0. Since for a small 
d transition matrices of the maps gd,ii, Qdj and gd,i9d,ii satisfy conditions of lemma [71 this 
bound implies ( ITS]) . If 02 > CI3 (the proof for this case is not presented), then 

\g''{w, zi, Z2)\ < 9dj9dji9%{w, zi, Z2) where ki + k2 + h = k 

for any (w, 21,2:2) G any < and a sufficiently small e{d) > 0. To prove f l76l) . we 

examine the location of the sets (/fij and guj. 

• We prove ( 175|) for (w, 21, 2:2) G w/,/. In ujij{e) for a sufficiently small e, g{w,zi,Z2) G w/,/ 
(see (162|) ) and the map (7 is bounded from above by gd,i for d < 0. Hence, by lemma [TO] 
(applicable by virtue of (1721) ), for a sufficiently small d the map g'' is bounded from above 
in coi,i{e) by (7^^. By ([721) . the map (77- (and hence (^^j for a sufficiently small d) satisfies the 
conditions of lemma [7] for asymptotic stability of {w,Zi,Z2) = 0. By this lemma applied to 
gdj, for any S > there exists e/ > such that 

for any {w, Zi, Z2) G u;/,/(e/), ^^^^ 
\g^{w,zi,Z2)\ < g^j{w,zi,Z2) < 6 for any k>0 and limfc^oo fi'''('W, 2^1, 2^2) = 0. 

• Now we prove f[75l) for any (ti?, zi, Z2) G w. 

Depending on the values of ai, 02 and 03, there are three possibilities for the location of 
^(w//(e)): 

(i) g{(jjii{e)) C OJij for some e > 0; 

(ii) g{(jjjj{e)) (f. Ujj for any e > and {guii{S)) fl W// 7^ for any 5 > (02 > 03 is required 
for this case); 

(iii) (yf(w//(e)) ^ w/j for any e > 0; there exists Sq > such that g{ijjii{6)) fl w// = for all 
S > 60 (02 < 0,3 is required for this case). 

Proofs for these cases are similar, although differ in details. Only a proof for case (iii) is 
presented here. 

By definition of sets uj and properties (b)-(f) of g, when e is sufficiently small we have 
in case (iii): 

g{un{e)) C W/, g{ujin{t)) C Vti, g{uiv{e)) C Vti, giuJv{e)) C Ujj, 

g{uiji{e)) Cuj\ uiji, g'^{uni{e)) Cuj\ ujiji, g'^{uiv{e)) C uij. 

For a small enough e, the map g is bounded from above by gdj in fi/(e) and by gdji in 
uj{e) \ uj. Hence, if {w,Zi,Z2) G uj{e) \ uj, then the iterates g^{w, Zi, Z2) are bounded from 
above: 

gd,ii{gd,igd,nY''~^^^'^iw, zi, Z2) for odd fc, 
{gd,igdjiY''^{w, zi, Z2) for even k 
till 5''^(w, Zi, 2:2) gets into w/^/ (if this happens). If (w, 2:1, 2:2) G a;/_//(e), then the bounds 



!/(«;, 2:1, 2:2)1 < 



1^^(^,2:1,2:2)1 < 



{gd,igd,iiY^ ^^^'^gd,i{w,Zi,Z2) for odd /c, 
gd,ii{gd,iigd,iY''~'^^^^gd,i{w, zi, Z2) for even k 
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hold true while g''{w, zi, Z2) ^ w/,/. Thus by lemmas [71 l9]and[T0l for e/ defined in (177]) there 
exists eji > such that 

\g'iw,z,,Z2)\<ej (78) 

for any {w,zi,Z2) G w(e//) \ w/,/ while g''{w, Zi, Z2) ^ uij- If g''{w,zi,Z2) ^ w/,/ for any 
> 0, then 

lim /(^i;,Zi,Z2) = 0. (79) 

• Finally, we prove f l75p for {w,zi,Z2) ^ 00. For a sufficiently small e > 0, the map g is 
bounded from above by gd,i in Qi{e) and by gf^^// in \ Qj. Hence, by (|6T]) . if (w, zi, Z2) G 
B f:\uj, then either 

|5f(w;, 2:1,^2) I < gd,i{w,zi,Z2) and g{w,zi,Z2) e u 

or 

|5(''(u7,zi,Z2)| < gd,n{uj,zi,Z2) while ^(''"^(w, 2:1, 2:2) ^ 
Thus, there exists e/// > such that for any {w, Zi, Z2) E -Bejjj 

\g\w, zi, Z2) I < en while g^'^{w, zi, ^2) ^ (80) 

If g^{w^ zi, Z2) ^ CO for any > 0, then 

lim g''{w,Zi,Z2) = 0. (81) 

Together bounds fl77j) . fl78|l and flHOj) imply that for any (ti?, Zi, 2:2) G 5^^^^^ 

|(y'*^(w, 2:1, 2:2) I < 5 for any k <0, 
and relations ([77]), ([79]) and dH]) imply 

lim g^{w,Zi,Z2) = 0. 

Part (a) is proven. 

(b.l) Let us consider the case 

a2 > 0, 03 > 0, ai + a2 < 1. 

(The last inequality is equivalent to < 1.) We will show that there exists S > such that 
for almost all (w, 2:1, 2:2) 

\g''{w, zi, Z2)\ > 6 for some /c > 0. (82) 

• We prove (l82l) for {w, zi, zq) G w/,/. For any d > and a sufficiently small e, the map g is 
bounded from below in 0Ji^i{e) by (/rf,/. Since gujijie) C w/,/ for a small enough e, by lemmas 
[7[ and [To] there exists 5/ > such that 

|(7^(w, 2:1, 2:2)1 > bi for any (w, 2:1, 2:2) G W/^/ and some > 0. (83) 
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• We prove now (182|) for almost any {w,Zi,Z2) € u. As in the proof of part (a), the three 
possibihties (i), (ii) and (iii), hsted above, should be considered. We again present a proof 
for case (iii) only, for which 02 < as- 
Consider the partition 

1^// = U_p<p<pQ(/i, /2, hp), 

where 

A ~ zr, f2 ~ zl^^'^-'"'^ and hp = 

Each set Q{fi, /2, hp) has the boundaries hp) and g(/2. Zip) (see (EH]) and ( pU]) ). By our 
choice of p and inequalities ( 150|) and (153|) . gq{f2, hp; ei) C a;/,/ for a sufficiently small ei. In 
case (iii), we have for a small enough €2 

gqifi, hp] €2) C uiji and g'^q{fi, hp; €2) C w//. 

Properties (b) and (c) of the map g imply 

2 ,^ /? Tx ? a ~ 02(1 +P + a2a) + 1 + aga ~ 

g q{f, h) = q{f, h) with / ~ 2:2 , a = , - , ^ , h ~ ^2, 

ai(l + P + 020;) 

for any q{f, h) C Uji such that / ~ Z2. By part (a) of lemmalHlthere exists Kp > such that 
g^^q{fi, h) = g(/i, h) for /i ~ Z2 and d < (ai — 03)"^. Hence, for a small €3 q^^{fi, h; €3) C 
cj/^/, and by virtue of fl^ there exists €4 > such that for any {w, Zi, Z2) G uJii{e4) at least 
one of the following inclusions hold: 

g^"{w, Zi, Z2) G oJij or g^{w, Z\, Z2) G a;/// for some k < Kjj. (84) 

Here Kjj = max_p<p<p Kp + 2. 
Now consider the partition 

^^iii = U_p<p<p(5(/i, /2, Zip), 

where 

/i = (-5^2(1 - \z2n) , f2 = [-^{l + \z2nj , hp = 4^^. (85) 

We have gQ{fi, /2, hp) = Q(/i, /2, /i), where 

/i(^2) ~ 72(^2) ~ 4, a = ((l+p)(ai-a3)+a2)/ai and - |z2|^(»i-'^^)/''i < /i< |z2|-('^i-«3)/-i; 
g'^Qifu /2, /ip) = Q(/i, A, /i), where 

/i(z2) ~ z^, f2{z2) ~ -22 ? = ■^(ai - 03)^7^0^"^ + a2ar^ and h ~ 2:2. 
(This can be checked by substituting the expressions for /i, /2 and hp fl85|) into fHTl) . ) 
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For 1 < k<K = Kjj/2 + 1, denote 

Hp,k = g-^'ig^'Qifi, /2, hp) n uju) n uju, (86) 

Hk = ^-p<p<pHp^k- (87) 
Note that Hp^k = Qifki, fk2, hk,p), where 

/fe, + CB'^Z2 ~ fk2 + 03-^x2 ~ N "fc > s, and /fc^ - fk^ ~ 7fc > "fc- 

By virtue of f lM|) . there exists €5 > such that for any (ty, 2:1, Z2) G ^///(es) at least one of 
the following inclusions hold true: 

g'^^{w, zi, Z2) G ojjj or g'^{w, zi, Z2) E uni for some k < 2K. (88) 

For a 6 > define the sets 

Hk{S) = HknBs, (89) 

Hk,M,-,kAS) = g-'''Hk2,...,kAS)nHk2,...,kM (90) 

and 

Gn{S) = Ui<fc (S). 

Let ee > be the e whose existence is asserted in lemma [5] for d = 1/2. By this lemma, for 
any 5 < eg 

Applying further lemma |5] (note that C, and Ck are independent of N and /Sm = /3i), 
we obtain ^ 



which implies 



Hence, 



Denote 



/^zi(G'7v((5)) < ^;U^,(GAr_i((5)). 



/i2i(G'7v(5)) < ^^izA^iii{6)) lim /i(Gjv(5)) = 0. 



Goo(^) = hm G^iS) 

N^oo 



(91) 



and set 6n = niini<j<6(ej). By definition of sets Gat (5), 

if {w, zi, Z2) e uiii \ GooiSii) and g''{w, zi, Z2) G Bsjj for all > 0, 
then there exists / > such that g''{w, Zi, Z2) G uij. 

Therefore, by (ES) and (ES]) 

for any (w, zi, ^2) G u; \ G'oo(57/7), \g^{w^ Zi, Z2)\ > 5iii for some k > 0, (92) 
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where Sni = min((57, Su). 

• Finally, we consider an arbitrary {w,zi,Z2). Let fiii{z2) and hiij{z2) be functions solving 
the system 

+ cf-^z2 = h{Dr + ETz^), Arh = f{Dr + etz^), 

which is equivalent to (7-invariance of q{fiii, hjji). There exists a small > such that 
(i) either q{fni, hjjj; 67) C fi/// or (ii) g(/m, hm; ej) fl Vlni = 0, 

depending on the values of ai, 02 and 03. 

Lemma [6](c) and property (d) of the map g imply that in case (i) there exists eg > such 
that for any {w, Zi, Z2) G VLjjj{t^) \ q{fin, hjjj) we can find k > such that (7^(x) ^ ^m- 

Define 

Qoo{S) = Uo<fc<oo5'"^G'oo(5) U qifiii, hjjj; 6ni) 

and note that /i(^oo(^m)) = 0. Set 60 = min(5///, 67, eg). By fl6T|) . fl62l) . fl92|) and definition 
of the set Qoo{S), 

for any (w, zi, 2:2) ^ ^oolf^o), Ig'^iu;, zi, Z2)\ > 5o for some A; > 0. (93) 

Since A* (^00(^0)) = 0, for 02 > and 03 > statement (b) is proven. 

(b.2) Next, we assume that 02 > and as < 0. By simple algebra (not presented), the 
latter inequality together with (l50l) imply < 1. 

• By the same arguments as in part (b.l), there exists 6j > such that for any {w, Zi, Z2) G 
ujij the statement fl83|) holds true. 

• Consider the partition 



^11— y^-p<p<pQifi, f2,hp) j U yJ-p<p<pQ{f3, fi,hp) 
where 

f — U / — U |-l/("i-a2) / _ 1^ |-l/(ai-a2) / _ 1^ 1-1/(01-03) U _ 

Jl — — P2I , /2 — — F2I , J3 — —\Z2\ /4 — — F2I ■,!>' — Z2 ■ 

There exists 5// > such that 

\g{w,zi,Z2)\ > Sii 

for any {w,Zi,Z2) G Q{fi, f2,hp,5jj). By part (c) of lemma [6], for any Q(f3, fi, hp) there 
exists f5{z2) such that 5'g(/5, /ip) = g(/5, /ip), and any {w, zi, Z2) E ujjj{ei) \ qif^, hp) satisfies 

\g''{w,zi,Z2)\ ^ uii 

for some k > 0. 

Define the sets Hk by (EE])- ([90]), 



,fcjv <oo -^^fci , . . . ,fcjv 

(6) (94) 
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and 

G^iS) = lim G^iS). (95) 

N^oo 

By the arguments employed in part (b.l), yu(G'oo) = and any {w,zi,Z2) G ujjj \ Goo has 
the property flOTl) : consequently, fl92l) holds true for some ^t-jt- > 0. 

• Define 

QooiS) = Uo<fe<oo^7~'Goo(5). (96) 

(now q{fiii, hjii) is not included into this set due to the bound 

\g{w,z,,Z2)\r^\{w,ZuZ2)r^^''^-'''^ 

for {w, zi, Z2) e qifiii, hjjj) and ai/ (ai - a^) < 1). For any {w, Zi, Z2) e ^00 ((^o) (where 
60 = min(5///,e8) and eg is defined in part (b.l) of the proof), the inequality \g\w, zi, Z2)\ > 60 
is satisfied for some / > 0. 

(b.3) Finally, we assume 02 < 0. Again, we study the sets guj and g^lj. The key points 
of the proof follow (details are omitted): 

• By ([HS]) and lemma [n](b), there exist 6jj > and Kjj > such that 

|5f''(x)| > 6jj for some < k < Ku, k>0 

for any x e cuij. 

• By property (b) of the map g and lemma [6](c) there exist a set 

q = U_p<p<pg(/p,o, hp) 
and e/ > such that for any {w, Zi, Z2) G W/(e/) \ g we can find > such that 

g''{w,Zi,Z2) i uJiiei). 



Consider the partition 
where 



UJjV — ^-p<p<pQ{fl, f2, hp), 



771 \ l/(ai-a3) / p \ l/(ai-a3) 

h={-j;^2{l-\z2m , f2=[-j;{l + \z2n) , hp = zl-'^. 



For 2 < A; < 00 define the sets 

Hp,k = g-^'ig^'Qifi, f2, hp) n ujv) n ujv (97) 

and the sets Hk, Hk{6), Hk, Gn{S) and G^oiS) by ([87D, (|89D, (190]), ^ and (|95]), respectively. 
There exists 6j > such that for any {w, Zi, Z2) & (Goo(^/) U q) 

\g''{w, z\, Z2)\ > 61 for some / > 0. 

• Define the set Qoo{S) by ( l96l) and note that for any {w, Zi, Z2) ^ (Gooi^i) U q) 

Ig^w, zi, Z2)\ > 6ni for some / > 0, 

where Sjjj < min(5/, 5//, e/). Since fi{Qoo{Si)) = 0, the theorem is proven. QED 
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Theorem 5 Consider the map 

g{w,Zi,Z2) = (Mkr',5u7"i +C^2k|"^L'w"^ + Ez2k|°^), (98) 
where ai > and 02(01 — 03) + (ai — 03)^ — ai > 0. Then 

(a) The fixed point {w, zi, Z2) = of the map g is asymptotically stable if and only if 

a2 > and > 0. 

(b) The fixed point {w,Zi,Z2) = of the map g is fragmentarily asymptotically stable if 
and only if 

a2 > —1 and 03 > 0. 

The proof is based on lemmas [7] and [H as the proof of theorem H] is based on lemma [71 
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